Helical glasses near ferromagnetic quantum criticality 
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We study the effects of quenched charge disorder on the phase reconstruction near itinerant ferro- 
magnetic quantum critical points in three spatial dimensions. Combining a replica disorder average 
with a fermionic version of the quantum order-by-disorder mechanism, we show that weak disorder 
destabilizes the ferromagnetic state and enhances the susceptibility towards incommensurate, spiral 
magnetic ordering. The Goldstone modes of the spiral phase are governed by a Sd-XY model. The 
induced disorder in the pitch of the spiral generates a random anisotropy for the Goldstone modes, 
inducing vortex lines in the phase of the helical order and rendering the magnetic correlations short 
ranged with a strongly anisotropic correlation length. 
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Quantum criticality remains one of the great enig- 
mas of modern condensed matter physics. Fluctua- 
tions around quantum critical points are now known to 
be responsible for many unexpected phenomena, such 
as the discontinuous quantum phase transitions seen in 
many itinerant ferromagnets (FM). Examples include 
MnSi [Hig, Sri_^Ca^Ru03 [2], C0O2 [3], UGes [4], and 
URhGe [5]. Such generic fluctuation-induced first-order 
behavior is a consequence of the coupling between the 
magnetic order-parameter- and soft electronic particle- 
hole fluctuations |6l [7] . 

Fluctuation-induced first-order behavior can be indica- 
tive of potential instabilities towards incommensurate or- 
dering [7]. Recently, this possibility has been explored 
within a fermionic version of the quantum order-by dis- 
order mechanism jSHTO]. The key idea is that certain 
deformations of the Fermi surface associated with the on- 
set of competing order enhance the phase space available 
for low-energy particle-hole fluctuations. The coupling 
of these additional soft modes to the order parameter 
fluctuations leads to a self-consistent stabilization of in- 
commensurate, spiral magnetic order, pre-empting the 
first-order transition to the homogenous FM. 

While disorder can smear continuous quantum phase 
transitions due to quantum Griffiths effects [11^ ^12, , the 
finite temperature tricritical point, below which the mag- 
netic transition turns first order and spiral order is pre- 
dicted, survives up to a critical disorder strength [6 . In 
this Letter, we argue that the interplay between disor- 
der and quantum fluctuations leads to a helical glass 
state consistent with the strongly inhomogeneous short- 
ranged ordered state observed recently near the avoided 
FM quantum critical point of CeFePO [13]. 

We consider electrons in three dimensions at chemical 
potential interacting via a local repulsion ^, and sub- 
ject to quenched charge disorder. The disorder potential 
Vcij) is Gaussian distributed with zero mean, variance 



and uncorrelated between different positions. The Hamil- 
tonian expressed in terms of the fermionic field operators 
^1){y) = (V^|(r),'0;(r)), is given by 

Our discussion of this model is divided into four parts. 
First, we summarize the phase diagram of the clean itin- 
erant FM. A fermionic quantum order- by-disorder treat- 
ment of fluctuations shows that FM order gives way to 
fluctuation- driven spiral order in the vicinity of the quan- 
tum critical point. Next, we treat the effects of weak 
disorder using the replica trick alongside the fermionic 
order-by-disorder approach. We find that although FM 
order is destabilized by disorder, the region of the phase 
diagram occupied by spiral order is enhanced. Next, we 
turn to the effects of disorder within the spiral phase. 
Taking account of lattice-induced anisotropy in the ori- 
entation of the spiral wave-vector and the combination 
of this with disorder, we show that the Goldstone modes 
are described by a 3d-XY model with random anisotropy. 
This form of disorder has a dramatic effect, driving the 
formation of vortex lines in the phase of the spiral and 
resulting in glassy behavior with short range magnetic 
correlations. Finally, we discuss our results in the con- 
text of recent experiments [13j and blue phases and lines 
of skyrmion defects predicted by other authors [Ml US] • 
Clean System: The free energy of the clean system de- 
scribed by ([1]) for Vc{r) = can, in mean field the- 
ory, be conveniently expressed in a Ginzburg-Landau 
expansion about the paramagnetic state: F^f{m) ~ 
am^ + bm^ + cm^ , with coefficients that are known func- 
tions of temperature and interaction strength. Restrict- 
ing for a moment to a spatially uniform FM, the fluctu- 
ation correction to the free energy takes the form jBUTO] 
Ff[{m) = \\g'^m^ In {nm^ + ^^), where A and n are con- 
stants [iF. Here and in the following we use rescaled. 
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FIG. 1: (color online) Fluctuations near the FM quantum crit- 
ical point stabilize an inhomogeneous, spiral magnetic state. 
While disorder (a^ = 0.4) destabilizes the FM, it enhances 
the region of spiral order. 



dimensionless units, T = T/fi and g = gpp {pf is the 
density of states at the Fermi level). This correction to 
the free energy shows a In T-divergence of the coeffi- 
cient, ultimately driving the transition first order. 

In fact, it has become apparent from subsequent analy- 
sis [SHTQ] that the first order FM transition is pre-empted 
by a transition into a spiral phase with order parameter 
m(r) = m[na^cos(q • r) + sin(q • r)] with q = qUz- 
For a spiral of pitch q (measured in units of the Fermi 
momentum) we may use the fact that the free energy 
is a functional of the mean-field electron dispersion - 
eiy(k) = k'^-\-h'^/ {kzqY + - together with the fact that 
the important integrals are peaked near the Fermi surface 
at low temperatures, to write the free energy as an angu- 
lar average over kz = kz/k^ F{m, q) = {F{^k1q^ + m?)). 

This results in relationships between powers of m and 
q in the Ginzburg-Landau expansion, e.g. the m^q^ co- 
efficient is proportional to the logarithmically divergent 
coefficient, explaining why the spiral is stabilized 
below some temperature. We expand in powers of q^ 
F{m, q) = /o(m^) + /2(^^)<7^ + ^f4{m'^)q^, but keep the 
full functional forms of the functions /2n(^^)- The phase 
diagram shown in Fig. [l] is obtained by minimizing the 
resulting free energy F{m^q). 

Replica Treatment of Weak Disorder: The effect of dis- 
order on the magnetic phase diagram, can be addressed 
using a combination of quantum order-by-disorder with 
the replica trick. We first recapitulate the steps in- 
volved in deriving the fluctuation corrections and incor- 
porate disorder at the appropriate stage. This proceeds 
in four steps: (i) decouple the interaction by a Hubbard- 
Stratonovich transformation in spin and charge channels. 



(ii) separate zero- and finite-frequency spin and charge 
fluctuations, allowing the possibility of a static magnetic 
spiral background, (iii) integrate out the electrons treat- 
ing disorder perturbatively to leading order. The latter 
introduces new terms in the free energy not found for the 
clean system, (iv) integrate out spin and charge fluctua- 
tions to quadratic order to obtain an effective Ginzburg- 
Landau description of the magnetic phase diagram. 

(i) The Hubbard- Stratonovich transformation is made 
by introducing classical fields p(r, r) and (/)(r,r), rep- 
resenting charge and spin, respectively. In order to 
compute the disorder averaged free energy we em- 
ploy the replica trick, F = —Tlimn^o{Z^ — l)/n. 
The effective replica action 5^, defined by Z'^ = 
J ?/^", p", 0"] exp(— S'n), consists of n copies of the 

clean system and a disorder vertex which is non-diagonal 
in the replica index and in imaginary time. 



5n = E/dr/ d\{g[{ry 



(2) 
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We have introduced the vector a = (cr^, cr^, a^)^ of Pauli 
matrices to conveniently express the electron spin. 

(ii) We decompose the fields p and (p into zero- and 
finite-frequency components, the former corresponding to 
static order in the system, p"(r,r) = po{r) + p"(r,r) 
and 4>^(y^t) = m(r) + 4>^(y^t). In the following, we do 
not consider the possibility of charge order (po = 0) and 
specify m(r) to be the spiral order parameter. 

Including m(r) in the free-fermion action allows a self- 
consistent free energy expansion in the presences of spiral 
magnetic order. Transforming the fermions to the rotat- 
ing spiral frame. 



(3) 



with cot6>(k) = kzq/m and Un denoting fermionic Mat- 
subara frequencies, we obtain a diagonal free-fermion 
action 5'o[(^, m] with corresponding Green's function 

(iii) Next, we integrate over the fermion fields (/?, 
treating the disorder vertex 6'dis[^,V^] perturbatively to 
leading order. The resulting mean-field action S'mf[m], 
which is independent of the fluctuation fields, is given by 



-(^dis[^,V^])o- (4) 



The average (. . .)o is taken with respect to the free- 
fermion action So[(p^(f^m] in the presence of spiral order. 
The fluctuation corrections, up to quadratic order in the 
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fluctuation flelds, are given by 



where {S'^S^w^) 



dis/ 



only includes 



connected diagrams. Here [p, 01 is the contribution 



arising from the flrst line of Eq. (l2|) while S[Lp^ Lp^ p, 0] 
denotes the linear coupling of the fluctuations to the 
fermion flelds originating from the second line of Eq. ^ . 

(iv) To obtain the fluctuation corrections to the dis- 
order averaged free energy, we perform the Gaussian in- 
tegrals over the fluctuation flelds, re-exponentiate, and 
flnally take the replica limit n ^ 0. 

Including disorder modifles both the mean fleld and 
fluctuation corrections to the free energy. The pro- 
portionality between expansion coefficients found in the 
clean system is broken in the presence of disorder since q 
enters not only the electron dispersion, but also through 
the transformation ([3| of the disorder vertex. The mean- 
fleld free energy contribution due to disorder is given by 



Tidis 



E 



(6) 



k,k' 



where Xuy{K^') = ^Eu; G^,{k,uJn)G^,^{k\uJn) = 
[nF{e,{k)) - nF{e,^{k'))]/[e,{k) - e,.{k')] with n^(e) = 

l/[exp (^^jr^ + 1] the Fermi- function and /i^y^-jy = 1 — 

hu,u ^ \{kz — k'^Yq^ /w?. The disorder strength in di- 
mensionless units is given by = cr^pi?//i. 

We have expanded hj^y up to quadratic order in q^ 
since we will only compute disorder corrections to the 
m^, m^, and q^m? coeflicients. Spin asymmetric combi- 
nations contribute to the latter, breaking the proportion- 
ality to the coeflicient. The leading disorder contri- 
butions to the fluctuation corrections are given by. 

Ft = r^lY. I Kyik,M){{-ir^'' 



( Xvy (ki , k2)x-u-u' (ks, ks 
- 2Kyy{ki,k2)n 

— v (k3)}, 



•k2) 



(7) 



with K,y{k^M) =fZ^^ G,(ki,a;,)G2,(k2,c.n) 

To obtain the corrections to the phase boundaries we 
compute and numerically and extract the con- 
tributions to the m^, m^, and g^m^ coefficients as func- 
tions of T and g. The resulting phase diagram is shown 
in Fig. [l] Disorder has two effects. On the one hand it 
destabilizes the homogeneous FM due to a positive con- 
tribution to the coefficient. On the other hand, it 
enhances the tendency towards spiral ordering because 
of a dominant, negative ffuctuation correction to the 
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FIG. 2: (color online) Three-dimensional plots of (a) the 
spiral ordering wave vector, (b) correlation length, and (c) 
anisotropy across the helical spin glass phase. 



q^rn^ term. While small disorder leads to a slight in- 
crease of the tricritical temperature Tc, this effect is re- 
versed by higher-order contributions that cut-off the non- 
analyticity of the clean system, Ffl ~ g'^m^ 

(/ the mean free path) and lead to a new non- 
analytic term of opposite sign [6l[T8]. Therefore, for suf- 
flciently strong disorder, the tricritical point is pushed to 
lower temperatures and is eventually destroyed. 



Disordered XY- Phenomenology in the Helical Spin Glass: 
Charge disorder does not destroy long range order (LRO) 
in the homogeneous FM. This is captured appropriately 
by our analysis of the replica averaged free energy. To 
see this, we may formally introduce independent order 
parameters in each replica. The resulting m^m| and 
(VmQ,)^(Vm/3)^ terms would arise from a replica average 
of disorder that couples to and (Vm)^, respectively. 
Such rotationally symmetric random mass and stiffness 
disorder does not destroy LRO in the FM. 

To understand the relevance of disorder in the ffuctu- 
ation driven spiral region we must analyze how disorder 
couples to the Goldstone modes. Taking into account lat- 
tice effects, the free energy is no longer invariant under 
continuous rotations of the spiral ordering wave vector 
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q. Starting, for example, from a tight-binding disper- 
sion on a cubic lattice and a chemical potential close to 
the bottom of the band, we obtain a small anisotropy 
term Fanis = A/ d^r[(5a,ma,)2 + (dyrriy)'^ + (d^m^)'^] 
(A > 0), which is minimized for q along one of the 
crystal axes. In the following we will assume q = quz 
as before. Note that the anisotropy term fixes the di- 
rection of q without changing the saddle-point equation 
= — /2(^^)//4(^^) for the modulus of q. The Gold- 
stone modes correspond to smooth deformations (j){r) 
of the phase of the spiral order parameter, m(r) = 
m[nx cos{qz -\-ny sm{qz + ^)], captured by the clas- 
sical anisotropic 3d-XY action 

\ i=x,y J 

(8) 

Disorder in the spin stiffness induces local random 
changes of q. While this obviously leads to disorder in 
the stiffness of the Goldstone modes it also generates a 
random anisotropy for the phase angle of the spiral, 

= ^Am2 I d\g{r) cos[20 + a(r)], (9) 

with a(r) ~ Sq{T)z a random phase and ^(r) = [{dya)'^ — 
(9cc<^)^]/4. No such random anisotropy is induced in the 
homogenous FM phase, leading to a very different re- 
sponse of the spiral phase to disorder. In particular, it is 
established that arbitrarily small disorder of random field 
or higher-order random anisotropy type destroys LRO in 
dimensions d < 4 [19H2T]. It is very likely that the ran- 
dom anisotropy 3d-XY model - and so the spiral phase - 
supports only short-ranged magnetic correlations [22] . 

This phase of short-ranged magnetic order is markedly 
anisotropic. In the following, we use the result by Dot- 
senko and Feigelman who argued that on scales much 
larger than the size of a vortex-line core, the correla- 
tion length of the isotropic 3d-XY model with random 
anisotropy disorder, D cos{n<p + ^i'^))^ spin stiffness 

2 2 

Ps is given by ^ = ^^§j^, where Kn{T) = (cos(n0))o 
is computed with the wave-like excitation in the absence 
of disorder. Specifying to n = 2 and taking into account 
a rescaling of the z-coordinate, we obtain 

2^5 (jot)' ('»' 

for the correlation length along the directions perpendic- 
ular to q where is the variance of ^'(r). The anisotropy 
is given by ^z/S.xy = \/4/4(m^)g'^/(Am^). Since from 
minimization of the free energy m and q are obtained 
as a function of temperature T and electron repulsion g 
we can evaluate ^ = £^{T^g) over the spiral region. In 
Fig. |2] the evolution of and ^z/^xy is shown. 



At the transition to the homogeneous FM, the magne- 
tization remains finite while g -> 0, leading to a diver- 
gence of ^. Note that in the presence of a magnetic easy 
axis anisotropy, this transition becomes weakly first order 
[24 . The correlation length changes significantly over the 
spiral region and is strongly anisotropic with > ^^y^ 
except very close to the transition to the FM. In the clean 
limit, cr^ ^ 0, LRO is recovered in the spiral phase. 

Discussions and Conclusions - By combining a replica 
disorder average with the fermionic quantum order-by- 
disorder approach [10] we have demonstrated that weak 
charge disorder destabilizes the FM phase but enhances 
the susceptibility towards incommensurate spiral order- 
ing below the tricritrical point. While disorder does not 
destroy LRO in the FM phase, the correlation length 
in the spiral region is finite even for infinitesimal disor- 
der. This is a consequence of the interplay of the (cubic) 
crystal anisotropy and disorder in the spin- stiffness which 
generates a relevant random anisotropy term in the ef- 
fective 3d-XY model for the Goldstone modes. 

A helical glass with short-range order is consistent with 
the peculiar incommensurate state recently found close 
to the avoided FM quantum critical point of CeFePO 
[T3] . The /iSR results indicate that the sample-averaged 
muon asymmetry function is dominated by strongly inho- 
mogeneous spin fluctuations and shows the characteris- 
tic time- field scaling expected from glassy spin dynamics. 
The anisotropy of magnetic fluctuations and the lack of 
evidence for FM cluster formation in field cooled magne- 
tization measurements rule out conventional spin-glass 
behavior [l^- Both observations find a natural explana- 
tion in terms of short range spiral order. 

The presence of oriented vortex-line defects in the he- 
lical glass should lead to a non-Fermi-liquid T^/^ tem- 
perature dependence of the electrical resistivity, as ob- 
served in the partially ordered phase of the helimagnet 
MnSi [i^[26]. From a symmetry point of view, helical 
magnetic phases are analogous to the ordered phases of 
cholesteric liquid crystals [14] and the partial magnetic 
order in MnSi is reminiscent of blue phases of cholester- 
ics [26|, where the line defects correspond to columnar 
skyrmion spin textures [15 . The vortex lines that we 
discuss here are rather different from skyrmion line de- 
fects. The latter correspond to vortices in the orientation 
of the spiral wave-vector, whereas the vortices predicted 
here are in the pitch of the spiral order for a fixed di- 
rection of wave vector. Nevertheless, some properties of 
line defects are quite generic. In Ref. [27], it was shown 
that columnar fluctuations of any type of line defects lead 
to a T^/^ contribution to the electrical resistivity. This 
mechanism should equally apply to the extended phase 
defects that we discuss here. 
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